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Abstract 

High-density fermion matter is meta-stable due to the anomalous non- 
conservation of baryon and lepton numbers in the electroweak theory. The 
meta-stable state decays by penetrating the sphaleron barrier separating topo- 
logically different vacua. The decay happens locally, and results in an anni- 
hilation of twelve fermions, accompanied by a production of gauge and Higgs 
bosons. We find numerically local bounce solutions determining the decay rate, 
which are classical paths in imaginary time, connecting two adjacent topolog- 
ical sectors. We also follow the real-time evolution of the bosonic fields after 
the tunneling and analyze the spectrum of the created bosons. 
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1. It is well known that transitions between topologically different sectors of 
the electroweak theory are accompanied by a change of the fermion number 
The height of the separating sphaleron barrier is about -E S phai ~ 10 TeV 
[|| and accordingly the tunneling rate is extraordinarily small; it is exponen- 
tially suppressed by the factor exp(— 2Si nst ) ~ 10 -153 with the instanton action 
<Sinst = 87r 2 /g r2 , g ~ 0.67. This prevents baryon and lepton number violation 
from being an observable phenomenon under ordinary conditions. 

This suppression is less strong, or the transition is even unsuppressed at all, 
if the temperature || |], ||, the particle energy || [7|, or the fermion density 
0- * s l ar S e - Hence, fermion number violating processes might have played a 
role in the earlier history of the Universe or could perhaps be realized in future 
supercolliders. 

In this paper we assume to have zero (or very low) temperature but a 
macroscopic amount of fermions of a very high density in thermal equilibrium. 
We describe them by the chemical potential fj,; since we have zero temperature 
this is the energy up to which the fermionic levels are occupied. In a sphaleron 
transition one level crosses the gap, all others are shifted such that after the 
transition we have the same spectrum again. But now one more level is occupied 
or depleted, depending on whether the levels went up or down, hence the created 
or annihilated fermion has the energy //. 

The fermion number Nf of each doublet changes with the Chern-Simons 
number Nqs of the bosonic fields as AiVf = ANqs- It means that the potential 
energy of the gauge and Higgs bosons is shifted ||, 1C|: 



V£ t = V pat + nN CS . (i) 

Some remarks are necessary here: First, we have fixed the zero-point of the 
energy to the trivial vacuum with Ncs = 0- Second, we may well neglect the 
change of the Fermi energy as due to the creation or annihilation of a few 
fermions. Third, AiVf = AiVcs actually applies to all fermion doublets, which 
also might have different chemical potentials. Therefore in eq. (|l]) we imply fi 
to be the sum of all these single chemical potentials. 

The former degenerate ground states of the system with integer iVcs now 
have energy [iNqs an d hence are metastable (Fig. 1). They can decay by 
tunneling through the sphaleron barrier to the adjacent vacuum with lower 
energy. If n exceeds a certain value // crit , the states become even unstable. 



In [10], arguments were given that the tunneling transition rate corresponds 
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Fig. 1: Schematic plot of the tunneling process between two topological 
sectors. 

to the baryon number violation rate at high accelerator energies. Thus the 
results of the present paper may also have some relevance for the high energy 
transition rate. 

The decay of metastable "false" vacua was discussed by Coleman already in 
1977 [pl[| . The transition rate per volume can be calculated in the semiclassical 
WKB-approximation, it is of the form 

T/V = Aexp(-B) . (2) 

Here B is twice the Euclidean action of the classical path that connects the 
decaying state and the escape point at the other side of the barrier and mini- 
mizes this action. It is a movement in the inverted potential; the system leaves 
the initial state, reaches the escape point with zero kinetic energy and would 
then move back the same way. Coleman called it "the bounce" . The prefactor 
A can be found from the small oscillation determinant about a bounce, with 
one negative mode removed. In the following we present the first numerical 
determination of the bounce in electroweak theory. 

2. We consider the minimal version of the standard electroweak theory with 
one Higgs doublet in the limit of vanishing Weinberg angle. In terms of dimen- 
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sionless rescaled quantities the bosonic part of the Lagrangian is 
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with the covariant derivative D, 
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\A a a T a , the field strength 



F 



2— fl- 

\F* u T a = i[D„,D v ] , F£ v = d^A a u - d v Al + e ahc A\Al, and the Higgs 
? ). v = mn /mw is the ratio of Higgs and gauge boson masses. 
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The fermions are not considered explicitly here, but by the modification ([l]) of 
the potential energy of the bosonic fields. We work entirely in temporal gauge, 
Aq = 0. The potential and kinetic energy and the Chern-Simons index are 
mw 
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Ncs is ° m y well-defined if the configuration space can be identified with the 
sphere S3, which requires the fields to be continuous at infinity. We fix them 
to the trivial vacuum there (A^ = 0, = (jj)). In general, we will look for 
fields that minimize the action and expect them to possess higher symmetries 
than completely arbitrary fields. Therefore we assume the fields to have the 
spherical symmetry of the sphaleron: 



4°(*,r) 
$(t,r) 



h [Oai ~ n a ni) h n a rii D(t, r) 



H(t,r) + iG(t,r)n-0](°) . 



(5) 



We assume the fields to be continuous and differentiable everywhere and to 
yield finite potential and kinetic energy. This poses some additional conditions 
on the radial functions at r = 0. 

For negative chemical potential fi, the false vacuum with Acs = can decay 
to the adjacent topological sector with iVcs = 1> where the ground state has 
energy < 0. The minimal energy barrier between the two sectors is given by 
the configurations of Akiba, Kikuchi, and Yanagida (AKY) Q which are static 
fields that minimize the potential energy Vp Qt for given values of Ncs (Fig. 1). 

The barrier vanishes completely if \i exceeds a critical value /i crit . From the 
condition that for \x > /i crit the potential Vp 0t gets a negative mode around the 
vacuum, one finds ||, 10 1 

16vr2 (rt 
Mont = — y~ m w ■ (6) 
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It is difficult to determine the bounce by solving numerically the equations of 
motion as an initial value problem. If we tried to find the propagation in time 
with a Runge-Kutta method or something similar, the system would always 
fall into some abyss of the potential —V^ ot . Instead, we rather prefer to find a 
stationary point of the Euclidean action 

S E = f° dt(T kin + V£ t ) (7) 

directly, without using the equations of motion, to is the time when the system 
reaches the escape point; since the lower bound is — oo, we can choose an 
arbitrary value for to, just as we fixed the origin for our spherical ansatz (|5|) 
somewhere. The corresponding translational invariance of the action does not 
influence the factor B in eq. (||), but contributes a space-time volume which 
allows to have a transition probability per time and volume. 

Se is a functional of the five functions A, B, D, H,G of eq. (||), depending 
on radial distance r and time t. Since the bounce has infinite extension both in 
space and in time, we introduce new variables x and u which cover only finite 
intervals. In practice, we use for example 

r(x) = X r arctan (| x) and t(u) = At arctan (| u) . (8) 

Using ansatz (||) and the substitution (g) we get 

Se = ^ P° du f 1 dx — L 2 f A 2 + B 2 + 2r 2 (H 2 + G 2 + \b 2 )) (9) 

2-7T J-l Jo W<P Y V J 

+ (<pA' + BD) 2 + [tpB' - AD) 2 + 2r 2 (^H' + \GD) 2 

+ 2r 2 (<pG' - \HD) 2 + -^{A 2 + B 2 - I) 2 + (H 2 + G 2 ){A 2 + B 2 + 1) 

+ 2A(G 2 -H 2 ) -4BGH+ \u 2 r 2 (H 2 + G 2 - l) 2 

+ 2p[D(A 2 + B 2 - 1) + ipBA' - ip(A - l)B' 

with 

, = * = = , w = «(«) = (£) , (io) 

and the dot and prime mean 4- and 4-, respectively. In order to find a station- 
ary point of S# numerically, we use a relaxation method which was discussed 
by Adler and Piran in great detail fll2"|| . The functional is put on a grid, let 
us say of size n u x n x , and the values of the five functions at the grid nodes 
are considered as 5n u n x independent variables. Of course there is no unique 
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way to discretize a functional, basically we followed the suggestions of 1 12 ] . We 
postpone the details to a more comprehensive publication. After putting some 
initial configuration on the grid, we sweep over it by changing the single vari- 
ables one after another. Each variable is changed according to the first iteration 
step of a Newtonian algorithm, which would find the stationary point of Se if 
it was considered as a function of the variable in question alone and all others 
were constants. This method is somewhat different from a procedure presented 



in [13] where a similar problem in the context of techni-baryons in the Skyrme 
model is solved. As mentioned above, the fields at x = 1 (r = oo) are fixed 
to the trivial vacuum, and some conditions must be obeyed at the origin. We 
take care that the initial configuration fulfills these properties and that we do 
not loose them during the sweeps. 

In principle, one has to perform some hundred or thousand sweeps until the 
situation is stable and the resulting configuration is a stationary point of Se- 
But there are some difficulties which must be handled carefully to get proper 
results: 

The main problem is that V^ ot is not bound, it can be negative. Due to 
energy conservation the bounce itself cannot have positive potential energy 
— Vp Qt at any time, its total energy E tot = T^ n — V^ ot is constant and zero. But 
in its vicinity, one can construct configurations which have positive potential 
— Vp Qt for some time and which give a lower action than the bounce. So actually 
the bounce is just a saddle point of Se, it is a minimum only if we restrict the 
potential energy — V^ ot to be non-positive. So we enforce this restriction, simply 
by rejecting all changes of fields rendering — VJf ot to be positive when we sweep. 
There are more sophisticated ways to take into account invariances like energy 
conservation p4[| , but in our case the simple remedy proved to be most effective. 
Besides, we use energy conservation to check if a given configuration is already 
close to the bounce and if the numerical accuracy (e.g. grid size) is sufficient. 

Another point is that the solution is not unique: We want to describe the 
bounce between the topological sectors with Nqs = and 1; but this still 
allows so-called "small" time-independent gauge transformations which do not 
change the Chern-Simons number. We choose the Nqs = vacuum to be the 
trivial one, which fixes the gauge completely. Practically, we have a boundary 
condition for the u = — 1 edge of the grid, but due to the large differences in the 
true time t near this boundary, the field variables are only weakly coupled there. 
We often observe that the fields tend to another vacuum near the boundary, 
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influenced by the main configuration at larger t. In order to have a smooth 
behavior, we sometimes apply a gauge transformation to the whole grid except 
the close vicinity of the boundary which already is in the trivial vacuum. 

Another degree of freedom is the choice of to, the time when the bounce 
reaches the turning point, i.e. the potential energy becomes zero again. For 
the initial configuration this point lies on the right end of the grid, during the 
sweeps we observe that it moves slowly to the left. For t > to the potential 
energy remains zero then (in the frame of our numerical accuracy). Actually, 
this shift to the left never stops completely, but finally the system moves as a 
whole, it just follows the zero- mode due to the free choice of to- 

A more detailed description of the numerical technique will be given in a 
subsequent paper. 

In Table 1 we give our results for the bounce action and the Chern-Simons 
number of the escape point for v = 1 (i.e. ran = mw) an d several values of 
p. The cases p = and 1 were not treated numerically, but clearly for p — > 1 
the action vanishes since the barrier disappears. For p = it is known that we 
have the instanton action. Actually, the finite vacuum expectation value of the 
Higgs field destroys the scaling invariance of the instanton which we have in 
the case of pure gauge fields, so the bounce is an instanton in the limit of size 
zero here ([Uj] . Only the introduction of the second massive parameter p allows 
the bounce to have finite size again in the other cases. 
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0.0 


1.00 


1.00 












0.2 


0.82 


0.85 


22.9 


1.9 


0.370 


0.020 


0.391 


0.4 


0.58 


0.68 


46.2 


5.9 


0.727 


0.060 


0.786 


0.6 


0.34 


0.51 


67.1 


9.6 


1.093 


0.102 


1.195 


0.8 


0.14 


0.30 


102.6 


16.8 


1.417 


0.180 


1.597 


1.0 


0.00 


0.00 













Tab. 1: Bounce action, Chern-Simons number at the escape point, 
and the number of gauge and Higgs bosons and their energies after 
the tunneling for several chemical potentials p = p/p C rit] the energies 
are given in units of rriw = ^W 1 ■ The mass of the Higgs boson 
is fixed to ran = myy ■ 



One might expect that Se could have many different stationary points, 
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in other words, there could be several bounce solutions with different escape 
points and values of the action. As a matter of fact, for our numerics we used 
at least two distinct initial configurations for each value of p, but always found 
the same solution. 

Fig. 2 shows the time dependence of the potential energy —V£ ot and the 
Chern-Simons number Nqs, with to fixed to zero. In Fig. 3 we plotted —V^ ot 
versus Acs and can compare this to the minimal energy barrier. Obviously, the 
kinetic energy terms force the system to take a path which deviates quite a bit 
from a possible path through the AKY configurations minimizing the potential 
energy only. As one may expect, the deviation grows with the height of the 
barrier. 

3. Next we want to examine the behavior of the system after the tunneling. 
To this end, we take the (numerically found) escape point configuration and 
let it propagate in the real Minkowskian time. We work on a two-dimensional 
grid again, but at this time we can solve the equations of motion using a kind 
of Runge-Kutta method. This was done in [jn| [lTj with a slightly disturbed 
sphaleron as initial configuration; we checked that our algorithm gives the same 
results in this case. 

We have to distinguish two different kinds of behavior. The system can 
either stay in the topological sector where it came to after the tunneling or 
move classically over the barrier towards the next ground state with even lower 
energy. For \p\ < 0.2, the system contains less energy than required to cross 
the next barrier, so that it will definitely stay in its topological sector. For 
values of \p\ between about 0.2 and 0.9 the system has enough total energy so 
that in priciple it could move over the barrier. The system, however, stays in 
its sector, because the energy is dissipated along the exitation modes in such a 
way that the system does not find a collective path towards the next vacuum. 
Finally for \p\ > 0.9 the system actually moves classically into the next sector, 
and subsequently also crosses the following barriers, i.e. it keeps moving into 
the direction of growing Acs • 

In the case when the system stays in its topological sector, we see that after 
some short time the potential and kinetic energies become almost constant. 
For the transition described above, from Acs = to 1, we observe very small 
oscillations of T^ in around |//|/2 and of VT ot around — |/i|/2, which is in accor- 
dance with the virial theorem. The energy is concentrated near a sphere which 
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expands with about 95% of the speed of light. iVcs becomes almost constant 
after a while, with small variations; but this constant is not the integer num- 
ber which classifies the topological sector. Nevertheless, the fermion number is 
changed by one for each doublet: The equation AiVf = ANqs is n °t valid here 
because in the Minkowskian setup the non-vacuum fields cannot be restricted 
to a finite region . 

Finally, we find that the time- averaged fields (^r^ la'" dt w ^h b — > oo) 
form a static vacuum configuration with V po t = and integer iVcs> and 

that at larger times the fields can be considered as small fluctuations around 
that vacuum. This allows to analyze the particle contents of the bosonic fields, 
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Fig. 3: The potential energy — V^ ot versus the Chern-Simons number 
Ncs for various values of the chemical potential p = /x/Vcrit- The solid 
lines correspond to the configurations which minimize the action (0, ||) 
(bounce trajectories), the dashed lines are the minimal energy paths 
( AKY-configurations) . 



we can perform a spectral decomposition in the basis of free particle modes 
[O]. We apply a gauge transformation to the fields which converts Af,^ to 
the trivial vacuum first, hence we work with fluctuations a, b, d, h, g around the 
latter one. The basis is found by solving the linearized equations of motion; 
basically we follow [16| here, but use V£ ot instead of Vp t- The free particle 
modes in our spherical ansatz are 



o*(t,r) = l(f3Ut)-Pl(t))rji(kr) 

h(t,r) = (p k (t)-^ k (t)-^ k (t))rj 2 (kr) 

+ (^(t) + (3 k (t)+(3 k (t))rj (kr) 

d k (t,r) = (-2f3 k (t) + f3 k (t) + (3 k {t)) j 2 (kr) 

+ (^(t) + (3 k (t)+(3 k (t))j (kr) 
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h k (t,r) = h k (t)j (kr 
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with 



g k (t,r) = -l$(t)kh(kr) 

$(t) = $ sin(cuit + ca) , h k (t) = h k sin(m + 7) , (12) 

uj 2 = 1 + k 2 , cj 2 j2 = 1 ± 2pk + k 2 , Vt 2 = u 2 + k 2 . (13) 

The j s are the spherical Bessel functions. For each time t and momentum 
k we can get the functions /?J*(i) and h k (t) from radial "Fourier" integrals; a 
comparison with (12) yields the coefficients (5 k and /i fc . 



In the linearized theory, the total energy is 

/•OO /"OO 

E = Ew + Eh = mw / dk rii(k) oji(k) + m\v / dA; njj(k) 0,(k) (14) 

j JO 

with the spectral densities 

Qtt2 _ 

(15) 



n (A;) = 


Qtt 2 

-^(fc 2 + l)^ (/3o fc ) 


ni(fc) = 


g 2 k 2 "1(A) 


n 2 (fc) = 




njj(/c) = 





Fig. 4 shows the spectral densities for the gauge and Higgs particles, njy = 
n o + n i + ra 2 and n#. The total particle numbers 

rco roc 

Ny/ = I dk nw(k) and Njj = dk nn{k) (16) 
Jo Jo 

and the energies E\y and Eh are listed in Table 1. The good agreement of 
Ew\-Eh with the energy fx gained from the fermions shows that the description 
in terms of the linearized theory is fully justified. 

4. In summary, our work falls into two main parts, both of numerical character. 
First, we found the classical bounce trajectory between two topologically differ- 
ent sectors as a stationary point of the discretized Euclidean action on a grid in 
time and the radial space-coordinate. A careful treatment is required since this 
is a saddle point rather than a minimum, and numerical problems due to the 
gauge freedom must be handled properly. Second, we examine the real-time 
behavior of the fields after the tunneling, which gives some hints about the 
possible bosonic signature of baryon number violating processes. Our ansatz 
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assumed a background of high fermionic density, but the results should also be 
interesting for the case of few particles with very high energy. 
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ation methods (0]), and thank P.Pobylitsa, M.Polyakov, and V.Petrov for 
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